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Abstract — We find the capacity region of linear finite-field 
deterministic networks with many sources and one destination. 
Nodes in the network are subject to interference and broadcast 
constraints, specified by the linear finite-field deterministic model. 
Each node can inject its own information as well as relay 
other nodes' information. We show that the capacity region 
coincides with the cut-set region. Also, for a specific case of 
correlated sources we provide necessary and sufficient conditions 
for the sources transmissibility. Given the "deterministic model" 
approximation for the corresponding Gaussian network model, 
our results may be relevant to wireless sensor networks where the 
sensing nodes multiplex the relayed data from the other nodes 
with their own data, and where the goal is to decode all data at 
a single "collector" node. 

I. Introduction 

Wireless sensor networks (WSN) consist of many nodes 
with sensing, computation, and communication capabilities, 
sharing a common wireless communication channel. In a 
typical WSN configuration, a large number of nodes measure 
possibly correlated data and transmit to a single collector node. 
This network problem is referred to as the "sensor reachback 
problem" in Q). In many applications, nodes are energy- 
limited and the physical distance between each sensing node 
and the common destination makes the transmission difficult 
or (energy wise) expensive. We investigate such a scenario 
where communicating nodes cooperate with each other and 
act as relays in order to transport their own data along with 
the data from the other sensing nodes. 

Wireless channels differ from their wired line counterpart in 
two fundamental aspects. On one hand, the wireless channel 
is a broadcast (shared) medium and the signal from any 
transmitter is received by potentially many receivers. This is 
called broadcast constraint. On the other hand, any receiver 
observes the superposition (linear combination) of signals 
from possibly many transmitters. This is called interference 
constraint. The simultaneous presence of these two constraints 
makes a general wireless network quite difficult to analyze. 
The multiuser Gaussian channel that models a relay network, 
unfortunately, has so far escaped a sharp general characteri- 
zation, even in the simplest case of a Gaussian relay network 
with a single source, single destination and a single relay Q . 
The capacities of Gaussian relay channel and certain discrete 
relay channels are evaluated in J3) and a lower bound to the 
capacity of general relay channel is presented. In [4|, capacity 
is determined for a Gaussian relay network when the number 
of relays is asymptotically large. 



In O, a simpler deterministic channel is proposed. While 
this channel model is significantly simpler to analyze, it is 
able to capture the key aspects of broadcast and interference 
constraints. For this model, referred to as the linear finite-field 
deterministic model, [5] determines the capacity for a general 
relay network with one source and one destination, as well as 
the multicast capacity with one source, multiple destinations 
and common information only. Our contribution in this paper 
builds heavily on the results and techniques of |5| and can be 
regarded almost as a trivial extension thereof. Nevertheless, 
to the best of our knowledge and somehow surprisingly, this 
simple extension has not been reported before. 

We consider the "sensor reachback problem" [ 1] for a linear 
finite-field deterministic network with arbitrary topology, a sin- 
gle destination node and independent information at the source 
nodes. We show that the capacity region for this network is 
given by the cut-set bound and takes on a very simple and 
appealing closed-form expression. Also, for a specific sources 
correlation model, we find necessary and sufficient conditions 
for the sources transmissibility. This result reminds closely 
Theorem 1 of [1], with the following main differences: on 
one hand, the result of fl] is more general since it applies 
to general correlated discrete sources observed at the sensor 
nodes and general noisy channels. On the other hand, our result 
applies to networks with broadcast and interference constraints 
while the result of HI requires "orthogonal" channels, i.e., with 
neither broadcast nor interference constraints. 

We expect that the achievability technique for the Gaussian 
(noisy) relay network proposed in [5] can be generalized to the 
case of multiple independent sources and a single destination 
as examined in this paper, so that a scheme that achieves a 
bounded and fixed gap to the capacity region in the Gaussian 
case can be found. Also, we believe that a fixed-gap rate- 
distortion achievable region can be found using independent 
quantization and Slepian-Wolf binning for the case of cor- 
related Gaussian sources with mean-squared distortion and 
Gaussian noisy channels, at least for some specific source 
correlation model (see O), especially matched to the discrete 
correlated source model considered here. This, however, seems 
to be a far more involved result since event in the standard 
case of Gaussian/quadratic separated lossy encoding (that 
corresponds to the case where the communication network 
reduces to a set of orthogonal links from the sensor nodes 
to the destination), a general fixed-gap characterization of the 
rate-distortion region is missing J6|. 

In this work we limit ourselves to the linear finite-field 
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deterministic model and we leave the fixed-gap achievability 
for the Gaussian case to future work. 

II. Review of the deterministic linear finite-field 

MODEL 

In this section we briefly review the deterministic channel 
model proposed in [5| and used in this work. The received 
signal at each node is a deterministic function of the trans- 
mitted signal. This model focuses on the signals interaction 
rather than on the channel noise. In a Gaussian (real) network, 
a single link from node i to node j with SNR srir^ j has 



capacity Cij = \ log(l + snr.ij) w log v /sn?7J. Therefore, 
approximately, mj — [log ^/snr^j] bits per channel use 
can be sent reliably. In [5 | (see also references therein), the 
Gaussian channel is replaced by a finite-field deterministic 
model that reflects the above behavior. Namely, the transmitted 
signal amplitude is represented through its binarjQ expansion 
X = J27Li B £ 2 ~ e where B e € F 2 . At the receiver, all the 
input bits such that y/s!Wij2~ i > 1 (i.e., received "above the 
noise level") are perfectly decoded, while all those such that 
y /sn? r i~j2~ l < 1 (i.e., received "below the noise level") are 
completely lost. It follows that only the most significant bits 
(MSBs) can be reliably decoded, such that the capacity of the 
deterministic channel is given exactly by riij and it is achieved 
by letting B\, . . . , B ni i.i.d. Bernoulli-1/2. 

A linear finite-field deterministic relay network is defined 
as a directed acyclic graph Q = {V,£} such that the received 
signal at any node j £ V is given by 
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i£V:(i,j)£S 



(1) 



where yj,x.i £ F|, sum and products are defined over the 
vector space ¥\, and where 



•• 
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is a "down-shift" matrix. Notice that n^j < q indicates the 
deterministic channel capacity for the link as described 
before. Without loss of generality, the integer q can be set 
equal to the maximum of all {rii j : £ £}■ The broadcast 
constraint is captured by the fact that the input x.j for each 
node i is common to all channels £ £. 

In the case of single source (denoted by s) single destination 
(denoted by d), Theorem 4.3 of (5'j yields the capacity of linear 
finite-field deterministic relay networks in the form 



C= min rank {Gg s c } 



(2) 



where is the set of cuts S C V, ft c = V— S such that s £ S 
and d £ S c , and where Gs.s c is the transfer matrix for the cut 
(S,S C ), formally defined as follows. Let Af(i) denote the set 
of nodes j for which £ £ (this is the "fan-out" of node 

'The generalization to p-ary expansion is trivial. Here we focus on the 
binary expansion as in Q. 



i) and let V(j) denote the set of nodes i for which £ £ 
(this is the "fan-in" of node j). The transfer matrix Gs^e is 
defined as the matrix of the linear transformation between the 
transmitted vectors (channel inputs) of nodes /3; n (5) and the 
received vectors (channel outputs) of nodes /3 ut(«S), where 
the inner and outer boundaries /3i n (<S) and /3 ut(<5>) of S are 
defined as |7): 



and 



An (5) = {t £ S : Af(i) n S c £ 0} 



/3out(S)={jGS c :7>(j)nS^0} 



In words: /3i n («S) is the set of nodes of S with a direct link 
to nodes in S c , and /3 ut(«S) is the set of nodes in S c with a 
direct link from nodes in S. 

Going through the proof of Theorem 4.3 in J5] we notice 
that the "down-shift" structure for the individual channels 
is irrelevant. In fact, this structure is useful in making the 
connection between the linear finite-field model and the cor- 
responding Gaussian case. As a matter of fact, if the channel 
matrices S q ~ ni - j in the above model are replaced by general 
matrices Sij £ ¥ q 2 xq , the result © still holds. 

III. Main result 

In a linear finite-field deterministic network defined as 
above, let V = {1, . . . , N, d}, where node d denotes the 
common destination and all other nodes {1,...,N} have 
independent information to send to node d. For any integer 
T = 1, 2, ... we let Wj = {1, ... , \2 TR >~\ } denote the mes- 
sage set of node i = 1, . . . , N. A (T, Hi, ... , Rn) code for the 
network is defined by a sequence of strictly causal encoding 
functions ff ] : W 4 x F^ t_1) -S- F|, for t = 1, . . . , T and 
i = 1, . . . , N, such that the transmitted signal of node i at (dis- 
crete) time t is given by Xj[i] = /^'(w^v^l], . . . , y*[i - 1]), 
and by a decoding function g : F^ 9 — > VVi X • • • Wat, such 
that the set of decoded messages is given by (uii, . . . , wn) = 
g(y d [l],...,y d [T]). 

The average probability of error for such code is defined as 
P n (e) = P((Wi, . . . , W N ) ^ (Wu . . . , W N ), where the ran- 
dom variables Wi are independent and uniformly distributed 
on the corresponding message sets Wi. The rate A^-tuple 
(Ri,...,Rn) is achievable if there exists a sequence of 
(T, Rt, . . . , R N )-codes with P n (e) -> as T -> oo. The 
capacity region C of the network is the closure of the set of 
all achievable rates. With these definitions, we have: 

Theorem 1: The capacity region C of a linear finite-field 
deterministic network (V,£) with independent information at 
the nodes {!,..., N} and a single destination d is given by 



E 

ies 



Rt < rank{G 5!lS c}, V5C {!,..., N}. 



(3) 



Proof: The converse of © follows directly from the 
general cut-set bound and by the fact that, for the linear 
deterministic network model, uniform i.i.d. inputs maximize 
all cut-set values at once 0, Q, 0. 

For the direct part, we build an augmented network by 
introducing a virtual source node and by expanding the 
channel output alphabet of each node i — {!,..., N}. Let 
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{no,i : i — 1, ...,N} be arbitrary non-negative integers. 
The channel output alphabet of node i in the augmented 
network is given by ¥ q 2 +no -\ The virtual source node has 
no = Y^,iLi n os input bits, partitioned into N disjoint sets 
Ui of cardinality no,i f° r i = l,---,Af, respectively, such 
that the bits of subset Ui are sent directly to node i and are 
received at the top no,i MSB positions of the expanded channel 
output alphabet. Fig. [TJ shows an example of such network 
augmentation for a "diamond" network [5|. 




Fig. 1. A diamond network with a source node 1, two relay nodes 2 and 3 
and a common destination d is augmented by adding node and virtual links 
to nodes 1, 2 and 3. 



After introducing the virtual source node, the augmented 
linear finite-field deterministic network belongs to the class 
studied in [5 | with the minor difference that the channel linear 
transformations are not necessarily limited to "down-shifts". 
Nevertheless, as we observed before, Theorem 4.3 of J5] still 
applies. Letting Rq denote the rate from the virtual source 
node to the destination node d, we have that all rates R 
satisfying 



Rq < min 

(fio,f2S)eAd 



rank {Gn 0i nc} 



(4) 



are achievable, where is the set of all cuts (£lo, °f ^ ae 
augmented network such that 6 flo and d <G 

For any such set £lo we have that £Iq = S U {0}, for some 
SC{1,.,., N}. Consequently, we have that fig = S c , where 
S, S c are subsets as defined in the statement of Theorem Q] 
Since the links from to any nodes i £ {1,...,N} are 
orthogonal by construction (not subject to any broadcast or 
interference constraint), we have that Go ,n^ has a block- 
diagonal form where a block is given by Gs,s c (the links of 
the original network, corresponding to the cut (Oq , f2g) via 
the correspondence f^o O S defined above) and other blocks, 
denoted by Goj for all j E S c , have rank noj, respectively. 
By construction, there is no direct link between and d so, 
without loss of generality, we can assume rio.d = 0. The 



general form for Ga ,n^ is 
Gs,s° 



G 











where we have indicated S c 
have 





Go. 



o 

Gq^sc 
. , i|<s e | }■ Therefore, we 



rank{Gsi ,sig} = rank{Gs, 5 <:} + ^ n j (5) 



In particular, the cut fi = {0} yields 



Rq<^2 no,. 



(6) 



By letting this inequality hold with equality, and by replacing 
this into all other inequalities, we obtain the set of inequalities 

^no^ranklGs,^}, V S C {1, . . . , N} (7) 



ies 



where we used the fact that X^jLi n o,j ~~ Sjes c n °>i = 
J2ies n <M- 

Consider now the ensemble of augmented networks for 
which there exist integers {n o i : i = 1, . . . , N} that satisfy 
(0. For such networks, the rate Rq — Ylj=i n o,j i s achievable 
(by [5|) and therefore the individual rates Ri = no^ are 
achievable by the argument above. Finally, the closure of the 
convex hull of all individual rate vectors R = (no,i , . . . , no,iv) 
of such networks is achievable by time-sharing. It is immediate 
to see that this convex hull is provided by the inequalities (OH 



IV. A SPECIFIC example: diamond network 

In this section we work out a simple example and provide 
an explicit achievability strategy. Consider the "diamond" 
network shown in Fig. [T] with nodes {1,2,3, d} and links 
of capacity ni, 2 , n-1,3, n-2,d and n 3: d- In this case, TheoremQ] 
yields the capacity region C given by 

Ri + R 2 + R3 < max{n 2 ,d, n 3 ^ d } (8) 

R1+R2 < n 2id + n lt3 (9) 

R1+R3 < »3,<j + ni,2 (10) 

Ri < max{n^ 2 ,ni,3} (H) 

R2 < n 24 (12) 

R3 < m 4 (13) 

Next, we provide simple coding strategies that achieve all 
relevant vertices of C. Any point R e C can be obtained 
by suitable time-sharing of the vertices-achieving strategies. 
There are 24 possible orderings of the individual link capaci- 
ties 711,2, n i,3, n 2,d and U3,d. Due to symmetry, the regions for 
the case n^^ > n 2,d will be the mirror image of the regions 

2 Indeed, the inequalities (5) represent the convex relaxation of the integer 
constraints (6). 
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for the case ri2.d > «3.d- Therefore, we shall consider only 
the cases where n-i.d > n3.d- 

The remaining 12 cases have to be discussed individually. 
For example, let's focus on the case nz,d < "1,2 < ^1,3 < 
ri2.d- An example of the network for the choice of the link 
capacities n$ t d = T"i,2 = 2, 77,1,3 = 3,n 2 ,rf = 4 is 
given in Fig. [2] Fig. [3] shows qualitatively the shape of the 
capacity region in the three possible sub-cases of the link- 
capacity ordering n^d < 711,2 < "1,3 < ri2,d'- case 1) for 
"1,2 + n 3,d < 77,1,3; case 2) for 711,2 + ^3,d > ^i,3> and case 
3) for ni 2 +«3,d > ^2,ci- In all cases, the achievability of the 
vertices B and C of the region of Fig. [3] is trivial, since these 
correspond to vertices of the multi-access channel with node 
2 and 3 as transmitters and node d as receiver. 

Case 1). Vertex A has coordinates (R\ = 77,1.2, R2 — 77,2. d — 

77.1.2 — «3.d, R3 = «3.d) and can be achieved by letting node 

1 send ni 2 to node 2. Node 2 decodes and forwards these 
bits after multiplexing its own n 2 .d — n i,2 — n 3.d > bits 
in the MSB positions, such that node 3 can send n 3 d bits 
without interference from node 2. Vertex D has coordinates 
(i?i = 711,2 + n 3:d ,R 2 = n 2 ,d - - ™3 : d,i?3 = 0) and 
can be achieved by letting node 1 send 711,2 +n^,d bits. These 
can be all decoded by node 3, then node 3 can forward the 
bottom (least-significant) 77,3, d bits of node 1 to node d. Node 

2 decodes the top (most-significant) 77,1, 2 bits from node 1, and 
forwards them after multiplexing its own bits. 

Case 2). Vertices A, D and E have coordinates (R\ = 
n-i.2, R2 = n 2 ,d - "1,2 - n 3 , d , R3 = 7i 3 , d ), (i?i = 774,3, R 2 = 
n 2 ,d ~ n lt3 ,R 3 = 0) and (R x = n lj3 ,i? 2 = n 2 , d - n xa - 
n 3,d, R3 = n i.2 + n 3.d ^ tin), respectively. Vertex A can 
be achieved in the same way as in Case 1). Vertex D can be 
achieved by letting node 1 send 1213 bits to node 3. Node 

3 decodes and forwards the bottom 713^. Since in this case 
«i,2 > 7ii,3 — n3,d, node 2 can decode the top 711,3 — «3,d bits 
of node 1, and forwards them to node d after multiplexing its 
own ri2, d — ni,3 bits, using its ri2. d — "3,d MSBs. Vertex E 
can be achieved by letting node 1 transmit 7113 bits, where the 
top 77,1.2 of which are received by node 2. Node 3 forwards 
the bottom 71 13 — 71 12 bits of node 1, and multiplex its own 
n 3.d + ?l i.2 — "1.3 bits. Node 2 forwards the top ti 12 bits 
from node 1, by multiplexing its own n 2 .d — "1,2 — n 3.d bits, 
transmitting over its ri2,d — ^ MSBs. 

Case 3). Vertices A, D and E have coordinates (Ri = 
ri2,d - n 3 , d ,R2 — 0, i? 3 = U3,d), (Ri = ni,3,R 2 = n 2 , d - 
"-1,3,-^3 = 0) and (Ri = 711,3, R2 = 0, R3 = n 2 A - "1,3), 
respectively. Vertex A can be achieved by letting node 1 send 
^24 — "3,d bits to node 2. Since 77,2. d — "3.d < "1.2 these can 
be decoded and forwarded to node d in the MSB positions. 
Node 3 simply sends 713^ bits to node d without interfering 
with node 2. Vertex D is achieved by letting node 1 send 

77.1.3 bits. The top 77,1,3 — n 3, d of these are decoded by node 2 
and forwarded together with 77,2. d — "1.3 own bits. The bottom 
77,3, d bits of node 1 are decoded and forwarded by node 3. 
Finally, vertex E is achieved by letting node 1 send 77,1,3 
bits. The bottom 77,3. d — ri2, d + 77,1,3 of these are forwarded 
by node 3, after multiplexing its own n%, d — "1,3 bits. Since 
77-2, d — 77.3, d < "-1,2, node 2 can decode the top ri2,d — 77.3, d 
bits from node 1 and forward them to node d using its MSB 



positions. Other cases follow similarly and the whole capacity 
region is achieved by decode and forward. 




Fig. 2. The configuration of the diamond network in the example (Case (1) 
in Fig [3) 




Fig. 3. The capacity region of the diamond network in the example. 

V. Transmissibility for correlated sources 

Consider the case of a sensor network where the nodes 
{1, . . . , N} observe samples from a spatially-correlated, i.i.d. 
in time, discrete vector source U = (J7i, . . . , Un) (see the 
source model in JTJ). The goal is to reproduce the source 
blocks u[l], . . . , u[T] at the common destination node d. If 
the source blocks can be recovered at the destination with 
vanishing probability of error as T — > 00, the vector source is 
said to be transmissible. In the case of a network of orthogonal 
links with capacities this problem was solved in [1] and 
yields the necessary and sufficient transmissibility conditiorH 

H(Us\U S c)< V «S C {1, . . . , N}. (14) 

3 The notation Us = {Ui : i £ S} is standard. 
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From the system design viewpoint, the above result yields 
the optimality of the "separation" approach consisting of the 
concatenation of Slepian-Wolf coding for the source with 
routing and single-user channel coding for the network (TJ. 

With the same assumptions and linear finite-field determin- 
istic network defined before, we consider a specific model for 
the vector source as defined in (6). Let Uq be a non-negative 
integer, and let V 6 Fj be a random vector of uniform 
i.i.d. bits. For all i = 1, . . . , N, let Ui C {1, . . . , n } and 
define Ui € F^ 1 ' as the restriction of V to the components 
{Vi : £ E Ui] of V. Then, the correlation model for the source 
(Ui, . . . , Un) is reduced to the following "common bits" case: 
sources Ui and Uj have common part {Vi : £ £ UiDUj} while 
the bits Vi in Ui —Uj and in Uj —Ui are mutually independent. 
It follows that H(Ui\Uj) = \Ui\ - \U t n U 3 \. 

This source model is somehow "matched" to a correlated 
source defined over the reals in the following intuitive sense. 
Consider N — 2 and let U\ and {/a denote the binary 
quantization indices resulting from quantizing two correlated 
random variables A\ € R and A 2 G R using "embedded" 
scalar uniform quantizers with n bits, such that their first m 
MSBs are identical and their last n — to least significant bits 
(LSBs) are mutually independent. If A\,A% are marginally 
uniform and symmetric, U\ and {/ 2 are exactly obtained by 
defining V as above, with no = 2n — m independent bits, and 
letting Ui include the to MSBs and the first set of n — to LBSs 
of V, and Ui include the same to MSBs and the second set 
of n — m LBSs of V. This model trivially generalizes to the 
case of N correlated sources and is related to the Gaussian 
sources with "tree" dependency considered in (6). For the 
source model defined above we have the following simple 
result: 

Theorem 2: The vector source U = (Ui,...,Un) is 
transmissible over the linear finite-field deterministic network 
(V,£) if and only if 

H(U s \U S c) <rank{G 5 , S c}, V S C {1, . . . , N}. (15) 

Proof: Again, we consider an augmented network with 
a single source node denoted by 0, with tiq output bits that 
we denote by V. As before, subsets Ui of cardinalities ng^ of 
these bits are sent to nodes i, respectively. However, differently 
from before we choose the subsets Ui to overlap in accordance 
with the vector source model. For the augmented network, 
the rate R from the virtual source to the destination d must 
satisfy (0). In particular, choosing Hq = {0} we get Rq < no- 
Generalizing the proof of Theorem[T]to the case of overlapping 
sets {Ui}, we find that for any cut (fioj ^o) °f tne augmented 
network such that Oo = S U {0} and fig = with S C 
{1, . . . , N} we have 

rank{Gn ,ng} = rank{Gs,sc} + rank{G ,sc} 

where Go.s= is the linear transformation between the inputs 
V and the (augmented) channel outputs of nodes j e S c . By 
construction, the matrix Go.s^ is formed by linear independent 
columns for all bits Vi with £ 6 {Jj £ gcUj. Therefore, 



Since V is uniform i.i.d., we have Rq = Uq = H(V) = 
H(U). Replacing these equalities into the set of inequalities 
© and using the chain rule of entropy H(XJ) = H(Us\U$c) + 
H(Usc) we obtain that the conditions (fl3T l are sufficient for 
transmissibility. On the other hand, if a source as defined in 
our model was transmissible, then the set of conditions (fT~5T > 
must hold, otherwise the rate Rq of the corresponding single- 
source single destination augmented network would violate 
©. Hence, necessity also holds. ■ 

VI. Conclusions 

In this work we have characterized the capacity region for 
a linear finite-field deterministic network with independent 
information at all nodes and a single destination node. In our 
setup, all nodes may relay information from other nodes as 
well as inject their own information into the network. This 
may serve as a simplified model for a large WSN where 
sensing nodes cooperate with each other to send the collective 
data towards a single collector node. For a specific model 
of discrete binary source correlation at the nodes, we have 
also found necessary and sufficient conditions for the source 
transmissibility. Albeit restrictive, this correlation model may 
be useful (e.g., see [6|) as a simple discrete "equivalent" (up to 
some bounded mean-square distortion penalty) for a spatially- 
correlated real sources whose components are observed and 
encoded separately at the network nodes. 

Motivated by these results, it is natural to investigate the 
performance of achievability schemes based on the techniques 
as in [5 1 (for independent information) and separated quan- 
tization and Slepian-Wolf binning (for lossy transmission of 
correlated sources) in order to achieve the capacity region 
or the distortion region of actual WSN, within a bounded 
performance gap. 
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